1 1 1 

Nonclassicality of quantum excitation of classical coherent fields in thermal 
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The nonclassicality of photon-added coherent fields in the thermal channel is investigated 
by exploring the volume of the negative part of the Wigner function which reduces with 
the dissipative time. The Wigner functions become positive when the decay time exceeds a 
' threshold value. For the case of the single photon-added coherent state, we derive the exact 

threshold values of decay time in the thermal channel. For arbitrary partial negative Wigner 
distribution function, a generic analytical relation between the mean photon number of heat 
, bath and the threshold value of decay time is presented. Finally, the possible application of 

SPACSs in quantum computation has been briefly discussed. 
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I. INTRODUCTION 



Nonclassical optical fields play a crucial role in understanding fundamentals of quantum physics and have many 
applications in quantum information processing Recently, the preparations of nonclassical non-gaussian optical 
, fields have attracted much attention. Usually, the nonclassicality manifests itself in specific properties of quantum 
— i ■ statisticSjSUch as the antibunching Q, sub-poissonian photon statistics [H, squeezing in one of the quadratures of 
£N| \ the field partial negative Wigner distribution etc.. The interaction between the nonclassical optical fields and 
. their surrounding environment causes the dissipation or dephasing [6(, which ordinarily deteriorates the degree of 
nonclassicality. The Wigner function is a quasi-probability distribution, which fully describes the state of a quantum 
system in phase space. The partial negativity of the Wigner function is indeed a good indication of the highly 
nonclassical character of the optical fields. Reconstruction of the Wigner distribution in experiments with quantum 
Qh! tomography @, H, Q demonstrates appearance of its negative values. This peculiarity can not be explained in the 
J-i ■ framework of the probability theory, hence it does not have any classical counterparts. Thus, to seek certain relations 
between the partial negativity of the Wigner distribution and the intensity of thermal noise may be very desirable 
for experimentally quantifying the variation of nonclassicality of optical fields in thermal channels. The volume of 
negative part of the Wigner function has been proposed as a useful tool for characterizing the nonclassical optical 
fields 0,10. 

The photon-added coherent states (PACSs) were introduced by Agarwal and Tara The single photon-added 
coherent state (SPACS) has been experimentally prepared by Zavatta et al. and its nonclassical properties have 
been detected by homodyne tomography technology 0, [H[ . Such a state represents the intermediate non-Gaussian 
state between quantum Fock state and classical coherent state (with well-defined amplitude and phase) jlq ]. For 
the SPACS, a quantum to classical transition has been explicitly demonstrated by ultrafast time-domain quantum 
homodyne tomography technique. Previous numerical study has indicated that the partial negativity of the Wigner 
function of photon-added coherent states in photon loss channel will completely disappear as the decay time exceeds 
a threshold value [T3|. However, the exact threshold value of the decay time has not been explicitly given. In this 
paper, we further investigate the nonclassicality of photon-added coherent states in the thermal channel by exploring 
the partial negative Wigner distribution. The volume of the negative part of the Wigner function is analyzed. The 
threshold values of the decay time corresponding to the transition of the Wigner distribution from partial negative 
to completely positive are derived. For SPACS in thermal channel, it is shown that the threshold value of the decay 
time is independent of the initial seed beam intensity. In addition, for any initial partial negative Wigner functions 
in thermal channel, we also obtain a general formula about the threshold value of decay time which is given by 

jt c (n) — In e ~' l +2 n 2n i where jt c (n) denotes the threshold value of the decay time for the case of thermal channel 
with mean thermal photon number n. From this generic expression, we can clarify how the thermal noise shortens 
the threshold value of the decay time compared with the case of photon- loss. 

This paper is organized as follows: In Sec. II, we briefly outline the basic contents of the photon-added coherent 
states. The explicit analytical expression of time evolution Wigner function of the SPACS in the thermal channel 
is derived. In Sec. Ill, the dynamical behaviors of the volume of the negative part of the Wigner distribution for 
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the SPACS and the two-photon added coherent states (TPACSs) is numerically calculated. It is shown that the 
volume of the negative part decreases with the decay time. The volume of the negative part of Wigner function of 
the TPACS more rapidly decreases than the one of the SPACS, though their Wigner functions become non-negative 
at the same threshold decay time. A generic relation between the threshold decay time of arbitrary partial negative 
Wigner functions in the thermal channel and the mean thermal photon number of the thermal reservoir is presented. 
Moreover, the possible application of SPACSs in quantum computation has been briefly discussed. In Sec. IV, there 
are some concluding remarks. 



II. WIGNER FUNCTIONS OF PHOTO-ADDED COHERENT STATES IN THERMAL ENVIRONMENT 

Let us first briefly recall the definition of the photon-added coherent states (PACSs) [HI]. The PACSs are defined by 
^/ N ^ j a \ a )i where |a) is the coherent state with the amplitude a and (a) is the creation operator (annihilation 

operator) of the optical mode. N(a,m) = m\L m (— |a| 2 ), where L m (x) is the mth-order Laguerre polynomial. When 
the PACS evolves in the thermal channel, the evolution of the density matrix can be described by [|| 

dp 7(71+ 1) f f f 
— = (Zapa' — a' ap — pa' a) 

-\- — {2a) pa ~ aa) p — paa*), (1) 

where 7 represents dissipative coefficient and n denotes the mean thermal photon number of the heat bath. When 
n = 0, the Eq.(l) reduces to the master equation describing the photon-loss channel. 

The presence of negativity in the Wigner function of the optical field is the indicator of nonclassicality. For an 
optical field in the statep, the Wigner function, the Fourier transformation of characteristics function [l8| of the state 
can be derived by [HI, [2(| 

W(f3) = -Tr[(O e - d )b{p)pbt(p)\, (2) 

where O e = X^^Lo |2n)(2n| and = Y^Lo + ^)(^ n + 1| are the even and odd parity operators respectively. In 
the thermal channel described by the master equation (1), the time evolution Wigner function satisfies the following 
Fokker-Planck equation [2l[ 

where q and p represent the real part and imaginary part of 0, respectively. The time evolution Wigner function can 
be derived as following: 



00 />oo r\ —ft p\ —^yi 

I H q . i, - / ) = <-x P ( -. / 1 / / W T (x, y)W( q ~ V ^Zl —, V ,0)dxdy (4) 
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where 

W T (x,y) ^—^-eM- 2{X ^f ] ) (5) 
7r(l + 2n) 1 + 2n 

is the Wigner function of the thermal state with mean photon number n. Substituting the initial Wigner function of 
the pure SPACS [3 

W s ( Q ,p, 0) = -2^i(|2 g + ^-aP) e „ 2|g+ip _ a | 2 
TTLi{~\a\ 2 ) 



and the Wigner function of the pure TPACS [12 



w? M) = ^y. 2 ?:f\ -^-°r (7) 

7rL 2 (-|a| 2 ) 
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into the Eq.(4), we can obtain the corresponding time evolution Wigner functions. For the case of SPACS, the 
analytical time evolution Wigner function can be obtained 

2e 7t [(£ - c 2 Rea) 2 + (C - c 2 Ima) 2 + c 4 - 1] exp[-2( M 2 + u 2 )/{l + c 2 )] 

7r(l + H 2 )(l + c 2 )3 
[(exp( 7 i)-l)(l + 2n)] 1 /2 j 

Re(a) — <7exp( 7 i/2), 
Im(a) — pexp(7</2), 
Re(a) - 2gexp(7t/2), 

Im(a) -2pexp(7t/2). (8) 

In Fig.l, the Wigner function of the SPACS with a = 0.5 in the thermal channel with n = 1 at two different values 
of decay time are plotted. It is shown that the phase space Wigner distribution of the pure SPACS with a = 0.5 
loses its circular symmetry and moves away from the origin due to the appearance of a definite phase. The partial 
negativity of the Wigner function indicates the nonclassical nature of the single quantum excitation of the classical 
coherent field. The thermal noise causes the disappearance of the partial negativity of the Wigner function if the 
decay time exceeds a threshold value. The tilted ringlike wings in the distribution gradually start to disappear and 
the distribution becomes more and more similar to the Gaussian typical of a thermal state. Eq.(8) also indicates 
that the negative region of the Wigner function of the slightly decayed SPACS is always a circle. Fig. 2 illustrates 
the phase space Wigner distributions of the TPACS with a = 0.5 in the thermal channel with n — 1. The upper 
figure of Fig. 2 shows the region of the negative part of the Wigner function of the pure TPACS is a ring belt which 
is different from the case of SPACSs. The lower figure of Fig. 2 shows that the absolute value of negative minimum 
of the Wigner distribution decreases as jt increases. When yt exceeds a threshold value, the partial negativity of the 
Wigner distribution of the TPACS also completely disappears. 
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III. INFLUENCE OF THERMAL NOISE ON THE NONCLASSICALITY OF QUANTUM EXCITATION 

OF CLASSICAL COHERENT FIELDS 

Recently, the volume Pnw of the negative part of the Wigner function has been suggested as a good choice for 
quantifying the nonclassicality (Tol . [ill . [l2j ■ Pnw is defined by 

Pnw = I / W(q,p)dqdp\, (9) 
J a 

where f2 is the negative Wigner distribution region. In Ref.jlT}, we have investigated Pnw of SPACS and TPACS in 
the photon-loss channel. It was shown that Pnw an d entanglement potential defined in Ref.[22| exhibit the consistent 
behaviors in short decay time. 

In this section, we bring our attention to the influence of thermal noise on the nonclassicality of the quantum 
excitation of classical coherent optical fields by calculating Pnw- In Eq.(8), we have obtained the analytical solution 
of the time evolving Wigner function of the SPACS in the thermal channel. Following similar lines, we can derive the 
one of the TPACS too. But the expressions are lengthy and do not exhibit a simple structure. 

Based on the analytical expressions of the Wigner functions, we have calculated the volume of the negativity in 
Wigner functions for both the SPACSs and the TPACSs. In Fig.3, P NW of SPACSs and TPACSs with the parameter 
a = 1.5 is plotted as the function of "ft for nine different values of the mean thermal photon number n of the thermal 
channel. It is shown that the thermal noise deteriorates the partial negativity, and Pnw monotonically decreases with 
the decay time. The larger the thermal noise, the more rapidly Pnw decreases, which implies that the nonclassicality 
of the optical fields are very fragile against the thermal noise. From Fig. 4, it can be observed that, the volume of the 
negative part of the TPACS's Wigner distribution is more fragile than the one of the SPACS's Wigner distribution 
against the thermal noise if the seed beam intensity \a\ 2 is not very large, though the initial pure TPACSs have larger 
values of Pnw than the pure SPACSs with the same value of a. 

The above results also indicate that Pnw becomes zero at a threshold decay time 7 i c which depends on the value 
of n. For the case of the SPACS in thermal channel, wc can directly derive the threshold decay time from the Eq.(8). 
7i c can be obtained as follows: 



(10) 
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which shows the threshold decay time is independent of the seed beam intensity \a\ 2 of the SPACSs. For the case of 
the TPACSs, the threshold decay times also satisfy the Eq.(lO). In Fig. 5, we have plotted both the exact analytical 
and the numerical results of the threshold decay time as the function of n for the case of SPACSs with a — 0.5. Well 
consistent between the analytical results and the numerical solutions is found. 

For arbitrary nonclassical optical fields which have the partial negative Wigner distribution function, there exists 
a relation between the mean photon number n of the thermal reservoir and the threshold decay time "ft c (n) beyond 
which their Wigner function become positive. 

, . , +2n 

where 7^(0) is the threshold decay time in the photon loss channel. In the derivation of the Eq.(ll), we have 
reformulated Eq.(4) as 

oo poo 



W{q,p,-yt) =exp( 7 t) / / W (x' , y')W(e^' /2 q' - \fe^' - lx' , e 7t ' ' /2 p' - \J7^' - ly',0)dx'dy', (12) 

where 
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W (x',y') = -exp(-2(a;' 2 + y' 2 )) 

X 



V 



y 

e lt/2 



y/l + (l + 2n)(eT'-l) 
P ~ v/l + (l + 2n)(e^-l) P 

7 t' = ln[l + (l + 2n)(e 7 *-l)]. (13) 

Obviously, the right side of the Eq.(12) represents a scaled time evolution Wigner function in the photon loss channel, 
i.e. n — 0. Therefore we have 

W^p^t)/^ 1 = iy(°) (<?>', 7 0/e 7t ', (14) 

where W^°\q,p, jt) is the time evolution Wigner function of the optical field in the photon loss channel. Based on 
Eq.(14) and the definition of Pnw m Eq.(9), we can prove 

PNw(jt) = P$l(jt?), (15) 

where Pj^^y^t') denotes the volume of the negative part of the photon-loss- induced time evolving Wigner function 
at the time t' of the optical field initially described by the wigner function W(q,p, 0). Assuming we have known 
the threshold decay time of arbitrary partial negative Wigner function in the photon loss channel, we can derive the 
corresponding threshold decay time of that Wigner function in the thermal channel via the Eq.(15). If P^w^') 



NW\ 

„7*c(0) 



+ 2n 



becomes zero at the decay time jt' = 7i c (0), it can be found that PNwilt) becomes zero at *yt = ln ; , ,, 
according to the relation jt' = ln[l + (1 + 2n)(e 7 ' — 1)]. Thus, we finish proving the relation in Eq.(ll) for arbitrary 
partial negative Wigner function. This relation may have some potential applications in the measurement of the 
temperature of the surrounding heat bath. In Ref.[17[, we have numerically investigated the entanglement potential 
and partial negativity of the Wigner function of SPACSs in photon-loss channel. The nonclassical properties of single- 
photon subtracted squeezed vacuum states in amplitude decay channel or in the phase diffusion channel have also 
been studied by Biswas and Agarwal [23] . The photon loss channel is only a specific example of the thermal channel 
at which the thermal reservoir is at zero-temperature. Present results not only generalize the results in Ref.[17[, and 
give out an analytical expression of the evolving Wigner function of SPACSs in thermal channel, but also find out the 
exact analytical relation between the mean thermal photon number of the thermal channel and the threshold time 
concerning the complete disappearance of partial negativity of the Wigner function of SPACSs. Furthermore, for 
any initial partial negative Wigner distribution function in thermal channel, we derive a general relation between the 
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threshold decay time and the mean thermal photon number of the thermal channel. Roughly speaking, the formula in 
Eq.(ll) implies that the higher the temperature of thermal channel, the more rapidly the nonclassicality of the optical 

fields with partial negative Wigner distribution decays. From Eq.(ll), it is also easy to know that ^t c (n) ~ ^ ^J -1 
when n > 1, and jt c (n) ~ 7* c (0) - 2(1 - e~ 7 * c(0) )n when n < 1. 

In the past five years, much attention has been focused on the application of coherent states in quantum computation 
|24l. l25l] . SPACSs can also find many applications in quantum information processes such as the quantum computation. 
Similar to the schemes in the literatures, one can regard the vacuum state 10) and the SPACS , a*\a) as two 

y/N(a,l) 

orthogonal states of an optical qubit |0}l and respectively. As an illustration, in this state encoding, a nontrivial 
two qubit gate can be implemented using only a beam splitter. Consider the beam splitter interaction described by 
the unitary transformation 

U BS =exp[#(a6 t + a t 6)], (16) 

where a and b represent the annihilation operators corresponding to two qubits |k) q and \x)b which take the states 
|0) or -j===aJ \a) . It is easy to derive that the output state produced by such an interaction is 

U BS \0)a\0)b = |0)„|0) 6 , 

Ubs— 7 ===a i \a) a \0) b = ] (cos(0)a t + ism(<j>)tf)\ cos(</>)a) a |z sm(0)a) b , 
y/N{a,l) y/N{a,l) 

Ubs / —^ =b' f \0) a \a) b = , = (cos(0)6 f + i sm{cj))a) )\i sm{<f))a) a \ cos{(f>)a) b , 
y/N(a,l) y/N(a,l) 

Ubs * 1 > o t 6 t |a) a |a) 6 = * - , (cos^a* + % sinM^XcosO^ + i sin(0)at)|e^ Q ) a |e^a) b , (17) 
The overlap between the output and input states can be expressed as follows: 

a<0|6<(]|C/ks|0)„|0)i = 1, 

a (a| fe (0|a(7 B sa t |a) a |0) fe 



N(a,l) 

-- ^-^ a (0\ b (a\bU B stf\0)a\a) b 
N(a, 1) 

COS 2 ((f>)\a\ 2 +COS0 laPfco^,!) 

N(a,l) 

1 

a (a\ b (a\abU B sa*tf\a) a \a) b 



N 2 {a,l 

= |a| 4 e 4 ^ + 2|a| 2 e 3 ^ + cos(2(/.) 2 | Q | 2(e **_ 1) 
N 2 {a,l) 6 



(18) 



If <fi is assumed to be sufficiently small such that </> 2 |a| 2 <C 1 but |a| is sufficiently large that <fi\a\ 2 is of order one. 
Then Eq.(18) can approximately become 

«<0| 6 (0|C^s|0>a|0) 6 = 1, 
1 a (a\ b (0\aU BS (J\a) a \0) b 



N(a,l 

■■ ^—ra(0\ b (a\bU B stf\0) a \a) b 
N(a,l) 

« 1, 

1 

1 6 (a | abU B sa* | a) a \ a) b 



N 2 {a,l) a ^ 

«e 2i<AH \ (19) 

Furthermore, the condition 2<t>\a\ 2 — it will guarantee the realization of a controlled phase gate by the beam splitter 
acted on the above encoding space. 

The results in this paper may be useful for investigating the influence of the thermal noise on the quantum infor- 
mation processes in which the photon added coherent state is used as a resource. It is also interesting to investigate 
the fidelity of the above controlled phase gates in thermal channel. The details will be discussed elsewhere. 
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IV. CONCLUSIONS 

In summary, we have investigated the non-classicality of photon excitation of classical coherent field in the thermal 
channel by exploring the partial negativity of the Wigner function. The total volume of the negative part defined by 
the absolute value of the integral of the Wigner function over the negative distribution region is calculated and the 
partial negativity of the Wigner function can not be observed when the decay time exceeds a threshold value which 
depends on the mean thermal photon number. For the cases in which the seed beam intensity |a| 2 is not very large, 
it is found that the TPACSs more rapidly lose their partial negativity of the Wigner distribution function than the 
SPACSs with the same seed beam intensity. For the case of SPACSs in thermal channel, the exact threshold value 
of the decay time beyond which the evolving Wigner function becomes positive is given as -ft c = ln( ^^ ). For the 
arbitrary nonclassical optical fields with partial negative Wigner function, we also present a generic relation between 
the threshold decay time of the thermal channel and the mean thermal photon number of the thermal reservoir under 
the assumption of the known threshold decay time of the photon-loss channel. Finally, the possible application of 
SPACSs in quantum computation has been briefly discussed. Present results in this paper may be useful for checking 
the influence of the thermal noise on the quantum information processes in which the photon added coherent state is 
used as a resource. 
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List of Figure Captions 

Fig.l. The Wigncr functions of the SPACS with a = 0.5 in thermal channel with n = 1 are depicted for two 
different values of decay time yt. 

Fig. 2. The Wigner functions of the TPACS with a = 0.5 in thermal channel with n = 1 arc depicted for two different 
values of decay time yt. 

Fig. 3. The Pnw of the SPACS and TPACS with a = 1.5 in thermal channel are depicted as the function of yt for 
different values of mean thermal photon number n. From top to bottom, n = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9. 
Fig. 4. Pnw is plotted as the function of the decay time yt for SPACS and TPACS with different values of a. (Solid 
line) from top to bottom, SPACS with a = 0.1,0.5,1.0,1.5, respectively; (Dash line) from top to bottom, TPACS 
with a = 0.1, 0.5, 1.0, 1.5, respectively, n — 0.5. 

Fig. 5. The threshold decay time yt c beyond which Pnw = is plotted as the function of the mean thermal photon 
number n for the case of the SPACS. (Solid square) numerical results; (Solid line) yt c = hiy±§^. It is shown that 
yt c ~ 0.5/n when «>1. The numerical calculations are based on the SPACS with a — 0.5. 




FIG. 1: The Wigner functions of the SPACS with a = 0.5 in thermal channel with n — 1 are depicted for two different values 
of decay time yt. 




FIG. 2: The Wigner functions of the TPACS with a = 0.5 in thermal channel with n = 1 are depicted for two different values 
of decay time yt. 




FIG. 3: The Pnw of the SPACS and TPACS with a — 1.5 in thermal channel are depicted as the function of 'yt for different 
values of mean thermal photon number n. From top to bottom, n = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9. 
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FIG. 4: Pnw is plotted as the function of the decay time yt for SPACS and TPACS with different values of a. (Solid line) from 
top to bottom, SPACS with a = O.f , 0.5, f .0, f .5, respectively; (Dash line) from top to bottom, TPACS with a = 0.1, 0.5, 1.0, 1.5, 
respectively, n = 0.5. 
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FIG. 5: The threshold decay time 'ytc beyond which Pnw — is plotted as the function of the mean thermal photon number 
n for the case of the SPACS. (Solid square) numerical results; (Solid line) jt c — In yqr§^. It is shown that jt c ~ 0.5/n when 
n >• 1. The numerical calculations are based on the SPACS with a = 0.5. 



